Deviation from Universality in Collisions of Ultracold 6 Li 2 Molecules 
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Collisions of 6 Li2 molecules reveal a striking deviation from universal predictions based on long- 
range van der Waals interactions. Li2 closed-channel molecules are formed in the highest vibrational 
state near a narrow Feshbach resonance, and decay via two-body collisions with Li2, Li, and Na. For 
Li2+Li2 and Li2+Na, the decay rates agree with the universal predictions of the quantum Langevin 
model. In contrast, the rate for Li2+Li is exceptionally small, with an upper bound ten times 
smaller than the universal prediction. This can be explained by the low density of available decay 
states in systems of light atoms [G. Quemener, J.-M. Launay, and P. Honvault, Phys. Rev. A 75, 
050701 (2007)], for which such collisions have not been studied before. 

PACS numbers: 03.75.-b, 34.20.Gj, 34.50.-s, 67.85.Lm 



Recent advances in the preparation of ultracold sam- 
ples of molecules are beginning to reveal how chemical 
reactions can occur in dramatically different ways at the 
quantum level compared to what happens in thermal en- 
sembles l]. Beyond seminal experiments demonstrating 
the effects of quantum statistics and induced dipole mo- 
ment on exothermic atom-exchange reactions of 40 K 87 Rb 
molecules J3, there is a range of new possibilities to be 
explored [l|, y-Q , such as controlling collisions through 
magnetic/electric field tuning to access scattering reso- 
nances or level crossings, understanding whether reac- 
tions can depend sharply on the specific quantum state 
of the collision partners or on the details of short-range 
inter-particle interactions, characterizing the outgoing 
states of reaction products, demonstrating coherent con- 
trol of reaction cross sections, and more. 

There exists a simple, universal description for two- 
body inelastic collisions and chemical reactions of an ul- 
tracold molecule with another molecule or atom 
This quantum Langevin model assumes a large number 
of available exit channels in the short-range part of the 
interaction potential, leading to a unit probability of loss 
there, and leaving the decay rate dependent on only the 
long-range van der Waals interaction between collision 
partners. It has been validated in various experimental 
settings JlJL, involving heavier alkali molecules like Rb2 
KRb 0, RbCs [13, and LiCs [H. 
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These universal collisions have a 100% probability of 
loss at short-range and therefore do not depend on details 
of the interaction potential there, such as scattering res- 
onances or reactivity determined by matrix elements be- 
tween quantum states. From a chemistry standpoint, it is 
thus more interesting to search for examples of collisions 
that deviate from universality. Such deviations should 
be more prominent in systems with low mass and conse- 
quently a low density of available decay states 14. lB|. 
making 6 Li2, consisting of the lightest alkali atoms, a 



uniquely suitable experimental system. We observe that 
two-body collisions of Li 2 in the highest vibrational state 
with free Li atoms deviates sharply from universality, as 
reflected in an exceptionally small two-body decay coeffi- 
cient. In contrast, the rates for both Li2+Li2 and Li 2 +Na 
collisions are universal. A recent experiment inferred the 
rate of Li 2 +Li decay from atomic three-body loss, but 
in a model-dependent way with uncertainty overlap pin g 
both our measurement and the universal prediction [16| . 

To our knowledge this is the first experimental real- 
ization of collisions with ultracold molecules where loss 
is described by physics beyond universal long-range van 
der Waals interactions [l7|. Earlier work by the Rice 
group reported a surprisingly low decay rate for Li2+Li2 
collisions [l8| . many orders of magnitude smaller than 
the universal prediction. Our measurement demonstrates 
that this rate is in fact universal, resolving a puzzle that 
has been prominent for the last decade. 

We sympathetically cool fermionic 6 Li with bosonic 
23 Na in a magnetic trap, as described in our earlier work 
[r|. The number balance as well as the final tempera- 
ture of the atoms can be adjusted by changing the Na 
evaporation endpoint and the initial loading times from 
the two atomic beams. We can completely evaporate 
Na, leaving a pure Li gas, or interrupt evaporation part- 
way to obtain a Na+Li mixture. At the end of magnetic 
trap evaporation, the atoms are transferred into a single- 
beam optical dipole trap with 5 W power and wave- 
length 1064 nm. Then we spin-flip Li and any remain- 
ing Na atoms (\F,m F ) = 3/2, 3/2) -> 1 1/2, 1/2) and 
2,2) — > |1, 1) respectively) with simultaneous Landau- 
Zener radio-frequency (rf) sweeps at 15 G. An equal su- 
perposition of the two lowest Li hyperfine states |1) and 
|2) (corresponding to 1 1/2, 1/2) and |l/2, — 1/2) at low 
field) is prepared by two non-adiabatic rf sweeps at 300 
G, separated by 10 ms. After holding for a further 500 
ms, the superposition becomes an incoherent mixture of 
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FIG. 1. Experimental sequence for molecule formation. Af- 
ter jumping near the B = 543 G resonance from above, we 
sweep the magnetic field B at 0.6 G/ms across a region 1.4 G 
wide. Immediately after the sweep, a short pulse of resonant 
imaging light removes free Li from the trap, leaving a pure 
gas of Li2. After a variable hold time, remaining molecules 
are detected by jumping back above resonance and imaging 
dissociated free atoms, or B can be kept below resonance to 
confirm that LI2 is invisible to the imaging light. 



|1) and |2). We then further evaporatively cool the Li, 
with any Na present, to T/Tp = 0.2, and transfer the 
atoms into a second, more weakly confining single-beam 
optical dipole trap parallel to the first one, with trap fre- 
quencies (u z ,v r ) — (21,480) Hz. Optimized formation of 
Li2 takes place in this second trap. 

Initial molecule formation experiments are done with 
Li |1) and |2) states only, without Na present. The ini- 
tial number in each state is 2 x 10 6 , corresponding to a 
peak density of 4 x 10 12 cm' 3 . A magnetic field sweep 
across the narrow 6 Li Feshbach resonance at B = 543 G 
[H| (Fig. [TJ) converts the atoms into diatomic 
molecules. These are closed-channel molecules in the 
highest vibrational state of the 1 £+ potential. 

As in our previous work [l9j . we unambiguously ob- 
serve the signature of molecule formation by applying a 
short blast of resonant imaging light to remove free Li in 
the |1) state from the trap. The light mass of Li means 
that it will be ejected from the trap after a single recoil, 
and moreover at 543 G the imaging transition is cycling, 
so a pulse duration of 20 fis is sufficient to leave no trace 
of 1 1) atoms. After the blast, imaging the |1) state while 
keeping B below resonance gives a negligible signal, since 
there Li2 is invisible to the imaging light. After switch- 
ing B above resonance to dissociate Li2, we image atoms 
in the |1) state as a measure of the molecule number, 
which for our optimized sweep parameters gives a forma- 
tion fraction of 10% or a molecule number of 2 x 10 5 . A 
second, independent confirmation of molecule formation 
is obtained by turning on a magnetic field gradient of 10 
G/cm for 6 ms while holding in-trap below resonance, 
which pushes free Li atoms away while leaving the spin- 
singlet molecules unaffected (Fig. [2]). 

For molecule decay measurements, we use two consec- 
utive blasts of imaging light, resonant with states |1) and 
|2) respectively, to remove free atoms in both hyperfine 
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FIG. 2. Images of Li2 molecules. Absorption images of (a) 
molecules dissociated into free atoms above the 543 G Fes- 
hbach resonance, (b) molecules held below resonance where 
they are invisible to the imaging light, and (c) separation of 
molecules (left) and atoms (right) in a magnetic field gradient. 



states from the trap immediately after the molecule for- 
mation sweep. This leaves a pure sample of Li2 molecules, 
which undergoes rapid initial decay from their vibra- 
tionally excited state, slowing down with increasing hold 
time in a way that is consistent with two-body decay 
from molecule- molecule collisions (Fig. [3]). The non- 
exponential nature of the decay rules out that lifetimes 
are limited by off-resonant excitations from the trapping 
laser, and we have also checked, by holding Li2 in trap 
at up to 8 G below resonance, that decay rates outside 
the coupling region around resonance are independent of 
magnetic field, as expected. We determine a two-body 
decay coefficient /?Li 2 +Li 2 = 

6(2) x lO" 10 cm 3 /s, with 
the uncertainty dominated by systematic errors in deter- 
mining densities. The quoted result comes from averag- 
ing multiple data sets, including those from a different 
crossed-beam trap geometry with higher initial densities. 

Applying the |1) and |2) state blast beams at the end of 
the hold time instead of immediately after the molecule 
formation sweep allows us to measure molecule decay in 
the presence of free Li atoms. The Li density is much 
higher than that of the molecules, so the presence of 
the atoms should significantly increase the decay rate. 
Surprisingly, we find that Li gives only a small, non- 
observable contribution to the decay [Fig. [Ha)] when 
compared to the decay from Fig. |3l corresponding to an 
upper bound /3Li 2 +Li < 5 x 10~ n cm 3 /s. There should 
be no Pauli suppression of collisions for closed-channel 
Li2 molecules with Li. We confirm this by checking that 
there is no enhancement of the decay after spin-flipping 
one component of Li from |2) to |3). 

When Na instead of Li atoms are trapped with Li2, 
significant enhancement of the loss rate is observed for 
similar initial atomic densities [Fig.[4jb)], corresponding 
to /3Li 2 + r N.a = 4(1) x 10~ 10 cm 3 /s. Na and Li only interact 
weakly [22j , thus the presence of Na has a negligible effect 
on Li2 molecule formation. The mixture is produced with 
the same temperature and initial Li2 density as in decay 
measurements done without Na. 
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FIG. 3. Lifetime of a pure sample of L12 molecules without 
free atoms. The solid line is an exponential fit up to 5 ms 
hold time, giving decay time constants of 8.7(5) ms, while the 
dashed line is a fit to a full two-body decay function. 

Two-body decay of molecules is described by 

= — /3Li 2 +Li 2 flLi, — /?Li,+Li™Li — /3ui 2 +Na?"lNa (1) 

"I.:.. - - - 

where n represents local densities. Experimentally, we 
measure the decay of total molecule number Nl; 2 rather 
than local densities, so Eq. (jTJ) can be written, assuming 
Gaussian density distributions for L12, Li and Na 

23/2 Li2_ _ _ j3 U2+u ^nu 2 —Pu 2 +unu — /^Lb+Na^Na (2) 
^ V Li 2 

with n denoting peak in-trap densities. The various two- 
body decay coefficients /3 can thus be extracted by fitting 
exponential decay rates at short hold times and normal- 
izing by the initial peak densities. The factor 2 3 / 2 ac- 
counts for the variation of density across the trap. The 
effect of deviations of density profiles from Gaussian is 
much smaller than the quoted uncertainties for /?. 

Full expressions for trapped ideal Bose(Fermi) gases 
in the local density approximation are used to calculate 
peak densities nNa(Li) = ±( ) 3 ^ 2 Li 3 / 2 (±z) , where m 
is the mass of the Na(Li) atom, T is the temperature of 
the gas, and Li„(z) is the n-th order Polylogarithm [23j ]. 
We determine T before the molecule formation sweep by 
fitting Li time-of-iiight expanded 2D column density pro- 
files with the fugacity z — as a free parameter, giving 
T/T F = 0.2 or T = 400 nK. After the molecule formation 
sweep, flu is lower by a factor of two compared to be- 
fore the sweep, despite a molecule conversion efficiency of 
only 10%, because many more atoms are associated into 
molecules that are lost via collisions with other atoms in 
the time it takes to complete the sweep. 

For the IA2 density, the simplest assumption is that the 
density distribution is proportional to the Li density pro- 
file before the sweep, meaning that nLi 2 can be estimated 
from the ratio of total numbers 

nu-2 = n u N u . 2 /N u (3) 
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FIG. 4. Lifetime of Li2 in the presence of free atoms, (a) 
Decay with Li (solid squares) is almost identical to the decay 
of the pure L12 gas (open circles), (b) Decay with a similar 
initial density of Na (solid squares) is significantly faster than 
the decay of the pure L12 gas (open circles). 



and likewise for the reduced nu after the sweep. This 
is valid if we neglect the complicated density-dependence 
of molecule formation efficiency [24j and assume that Li 
atoms are selected at random from the Fermi sea to form 
molecules that do not have time to reach thermal equi- 
librium. We improve on this by accounting for the ef- 
fect of equilibration. Molecules form with the same aver- 
age center-of-mass kinetic energy as the free atoms, but 
twice the potential energy (due to their larger polarizabil- 
ity). Assuming that equilibration distributes this excess 
energy among all the degrees of freedom (according to 
the Virial theorem applied to harmonic traps) , the cloud 
width is rescaled by a/3/4. In our experiment, L12 does 
not equilibrate along the weak axial trapping direction, so 
the radial cloud diameter is instead rescaled by a/7/10. 
This implies that Eq. §3§ underestimates the peak L12 
density by about 30%, which is within the quoted un- 
certainty of our results, but we nevertheless include the 
correction in the analysis. This correction is unnecessary 
for hu after the sweep [25|. 

Figure [5] shows how (3u 2 +u deviates sharply from from 
the universal prediction while /3u 2 +u 2 an d /?Li 2 +Na are 
universal. In the quantum Langevin model for universal 
collisions of ultracold molecules, the assumption of total 
loss at short-range leaves the decay rate dependent only 
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FIG. 5. Predictions for two-body inelastic decay coefficients 
from the universal model (open circles) compared to experi- 
mental measurements (vertical bars). 



on the long-range van der Waals interaction via 0, IH 



P = 9- 



fi 



0.48( 



(4) 



where the prefactor g = 1(2) for (in) distinguishable colli- 
sion partners accounts for wave- function symmetrization, 
fi is the reduced mass of the two-body system, and with a 
length scale a for van der Waals interactions expressed in 
terms of the Cq coefficient. The weak C^ 4 dependence 
means that variations of /3 come mainly from differences 
in \i. For collisions involving weakly bound molecules, Cq 
can be approximated as the sum of the correspondin g co - 
efficients for all combinations of atom pairs involved 13] , 
with the atomic Cq coefficients taken from calculations 

The observed deviation from universality in Li2+Li 
collisions can be explained by a quantum depletion 
model, which predicts a two-body coefficient of 5 x 10 -11 
cm 3 /s for Li 2 in the least bound triplet 3 S+ state [28j |. 
and a similarly small value in the low 10 -11 cm 3 /s range 
for the least bound singlet state [14[ that is con- 
sistent with the measured value of f3u 2 +u- This effect 
is absent for Li2+Na, Li2+Li2, and collisions involving 
heavier alkali molecules, because these all have a higher 
density of decay states in the exit channel [l4 , 28 1 . The 



effect is also only present for the highest vibrational state, 
in which the atoms spend the longest time near the outer 
turning point of the van der Waals potential, where they 
do not have sufficient kinetic energy to transfer to a col- 
lision partner in a vibrational relaxation process. Li 2 in 
the second highest vibrational state is exp ected to have 
an order of magnitude larger f3u 2 +u [28(. Confirming 
this prediction requires using a two-photon Raman tran- 
sition to change the vibrational quantum number. 

Closed channel molecules formed around a Feshbach 
resonance in the highest vibrational state have size on 
the order of the van der Waals length a 0, [13], while 
the universal model assumes a clean separation between 
loss processes at short range and the long-range van der 
Waals interaction. Our measurements show that Li2+Li2 



and Li2+Na collisions are still well described by universal 
redictions. This follows previous work, both theoretical 
and experimental jlQj, indicating that decay rates are 
mostly independent of the vibrational quantum number 
of the molecule. Decay of other molecules like Na2 [HI, 
30l ] and Cs2 31 1 from their highest vibrational states also 
show fair agreement with universal predictions. 

Our work resolves the puzzle of unexpectedly long life- 
times of Li 2 molecules [18J observed by the Rice group 
[32l |. This was one of two seminal reports of long-lived 
fermion pairs H, 33| that triggered the subsequent explo- 

While the obser- 
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ration of the BEC-BCS crossover 
vation of long lifetimes near a broad Feshbach resonance 
[33| was quickly explained in terms of Pauli suppression 
of collisions involving open-channel-dominated fermion 
the other observation near a narrow reso- 
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pairs 

nance has remained unexplained, since Pauli suppression 
should be absent for closed-channel molecules outside 
the narrow coupling region around resonance. Our mea- 
surements demonstrate conclusively that the lifetimes of 
these molecules are in fact short, as expected (Fig. [3]). 
The atomic densities and temperatures in our experiment 
are comparable to those in the previous report. In ad- 
dition, we have reproduced their magnetic field sweep 
during the hold time [35| instead of the usual holding 
at a constant field, and found no significant lifetime en- 
hancement. Note that the Rice experiment deduced the 
presence of molecules only from differences in Li atom 
numbers [l8j . whereas our experiment unambiguously 
identifies Li2 molecules through magnetic field gradient 
separation, and by their survival when Li is removed by 
resonant imaging light (Fig. [5]). The Rice results are 
also incompatible with recent work inferring /?Li 2 +Li from 
three-body atomic loss [lfj} , but this has not been pointed 
out before. 

In summary, we have observed that, for Li2 molecules 
in the highest vibrational state formed around a nar- 
row Feshbach resonnce, Li2+Li collisions deviate sharply 
from the universal predictions of the simple quantum 
Langevin model, while Li 2 +Li2 and Li2+Na collisions are 
universal. This is the first example of collisions involv- 
ing ultracold molecules with loss determined by physics 
beyond long-range van der Waals interactions. 
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